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FREE -STREAM TEMPERATURE, DENSITY, AND 
PRESSURE MEASUREMENTS IN AN 
EXPANSION TUBE FLOW 

By Kenneth V. Haggard 
Langley Research Center 

SUMMARY 

An experimental study has been conducted to determine test-flow conditions in the 
Langley pilot model expansion tube. Measurements of temperature, density, wall pres- 
sure, pitot pressure, and shock and interface velocities have been compared with theoret- 
ical calculations based on various models of the flow cycle. The vibrational temperature 
and integrated density of the molecular oxygen component of the flow were measured by 
use of vacuum ultraviolet absorption techniques. These measurements indicate both the 
presence and possible degree of nonequilibrium in the flow. Data are compared with 
several simplified models of the flow cycle, and data trends are discussed. 

INTRODUCTION 

The Langley pilot model expansion tube is a ground based facility capable of pro- 

f 

ducing a short-duration flow with a velocity, density, and temperature which simulate the 
conditions encountered by an entry vehicle. The method used to produce the test flow not 
only can be used with a variety of test gases but also avoids the creation of the high- 
enthalpy stagnation region present in most other entry simulation devices. The elimina- 
tion of the stagnation region reduces the effects of chemical and thermal nonequilibrium 
on the flow cycle. 

Previous experimental studies of expansion tube flows (refs. 1 and 2) have shown 
that the wall pressure and pitot pressure in the test flow are significantly different than 
those theoretically predicted. Although the importance of nonequilibrium effects is 
reduced in a facility of this type, its effects can alter the test condition. Thus, the pres- 
ent study was undertaken to characterize the degree of nonequilibrium in. the flow by meas 
uring the vibrational temperature of the test flow in the Langley pilot model expansion 
tube. The experimental method used the well-known vacuum ultraviolet absorption prop- _ 
erties of the oxygen molecule and constituted a measurement of the vibrational tempera- 
ture of the oxygen molecule. Integrated oxygen densities along the absorption path are 


V 


velocity, m/s 


velocity of secondary interface, m/s 
Vg i velocity of primary shock, m/s 

x distance along tube, m 

Z compressibility factor, p/pRT 

a = 147 (see eq. (4)) 

A 160 

y lsentropic exponent 

p tc 

A = -p 2 - (see eq. (9)) 

X wavelength, nm 

p density, kg/m3 

p J scale depth, cm (see eq. (6)) 

d 

p density of C>2, kg/m 3 

°2 

Subscripts: 

0 standard conditions 

1 initial state in intermediate chamber 

2 state behind primary shock 

5 state of test- gas flow 

10 initial state in acceleration chamber 
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These flow models neglect the effects of diaphragm losses, viscosity, extraneous 
shock waves, and finite chemical and vibrational relaxation rates. References 3,4, and 
5 give calculated performance curves based on these models. The effects of finite chemi- 
cal and vibrational relaxation rates have been considered in references 6, 7, and 8. The 
flow model used in these reports was inviscid, free of reflected shocks, and assumed equi- 
librium behind the primary shock, while the chemical model contained only N 2 , O 2 , N, 
and O. These restrictions and the uncertainty in the values of the chemical and vibra- 
tional relaxation rates limit the use of these calculations to providing a qualitative feel 
for the processes being studied. These nonequilibrium models, though simplified, do 
suggest that for the temperature and density ranges of concern in expansion- tube flows, 
the chemistry will depart from equilibrium in the test flow because of the rapid cooling 
in the expansion fan. 

Calculation of the test-flow conditions has frequently been based on the flow cycle 
of figure 1(a). For a primary shock wave of a given strength, the state in region 2 is 
determined from the Rankine-Hugoniot relations on the assumption of local thermody- 
namic equilibrium. An unsteady expansion at constant entropy then yields the velocity 
of the test flow as a function of pressure or any thermodynamic state property. 

It would be desirable for comparison with the present results if finite relaxation 
rates were included in the calculations so that nonequilibrium effects could be predicted. 
The only nonequilibrium calculation method available is that of references 6, 7, and 8. 
However, because of the limited utility of this program and the necessity to perform a 
separate calculation for each streamline, no finite-rate calculations were performed in 
the present study. Instead, a limiting situation of zero- relaxation rate (frozen flow) was 
calculated, in which the composition and vibrational energy are assumed to remain con- 
stant through the expansion at the equilibrium condition established behind the primary 
shock wave. This calculation takes the form of a "perfect-gas" isentropic expansion of 
a mixture of atoms (y = 5/3) and diatomic molecules (y = 7/5). The isentropic exponent 
of the mixture may be computed from the compressibility factor in state 2 behind the pri- 
mary shock by the relation 


4 + 3Z 2 

r 2 = T7z7 


(i) 


A different calculation method is to incorporate the effects of the reflected shock at 
the secondary diaphragm location. However, the flow situation sketched in figure 1(b) 
shows that the reflected shock varies in strength for the different stream paths, depend- ’ 
ing on the extent to which it has been weakened by interaction with the expansion fan, 
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The stations at which p w and p t were measured were located just upstream and 
just downstream of the ultraviolet beam. Time correlations between these various meas- 
urements were made by using the measured interface velocity. 

DATA REDUCTION 

The vibrational temperature and density of O 2 were determined from the absorption 
measurements by using the method of references 11 and 12. Basic data and a theoretical 
treatment are found in references 13 and 14. With this technique the absorptivity of the 
flow is measured at two wavelengths. These wavelengths are selected so that the absorp- 
tion coefficient at one is very sensitive to temperature whereas the absorption coefficient 
at the other is relatively insensitive to temperature. The degree of absorptivity is also 
controlled by the density of O 2 in the flow. Thus, since the absorption coefficients are 
functions of only temperature and wavelength, the temperature is easily determined. 

This procedure will become clear by considering the equations governing the absorption 
process. 

The absorption of the ultraviolet radiation was assumed to follow Beer's law. Thus, 


= h,o^ ex p 


j 5- J P(t)k x (t) d l 


( 2 ) 


where p is the absorber density and the integral is taken over the absorption path. The 
parameter A^ at A is 


A x (t) = log 


A.o 


(t) 


I>(t) 


rpWMt) 

J P n 


(3) 


and the ratio of A 147 and A-^qq is defined by 


g(t) ,^. W^ (4) 

A 160 (t) J p(t)k 16Q (t) d 1 


where subscripts 147 and 160 represent values of A. 


The data-reduction model for density and 
assumptions: 


v >°2 


includes the following 
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The result of this procedure is shown by the band in figure 5. The width of the band indi- 
cates the probable error in the average. The results of applying this averaging process 
to a series of different nominal conditions are shown in figure 6. 

Aside from the initial acceleration of the interface, there is no clear trend to indi- 
cate any systematic change of the velocity as the interface moves down the tube. The 
lack of a pronounced deceleration is in disagreement with the results of reference 1. 
These previous results are shown in figure 7. These data show a definite overshoot 
occurring 1.2 meters downstream of the secondary diaphragm, and the severity of the 
overshoot increases with increasing p^. An examination of the equipment and data- 
reduction procedure has not shown any reason for this discrepancy. The variation of 
the tube diameter was checked to determine whether significant variations in the length 
of the standing wave could occur, and it was found that they could not. Even if such a 
pattern had been found, it would not account for the observed dependence on p . Dif- 
ferent tube sections were used for these tests, but again it is not understood how this 
could produce such an effect. In view of the uniform measured during this study, 
no attempt has been made to use the interface attenuation model described in reference 1. 

The measurement of was standardized by using the value obtained from the 
microwave technique over the last 1.6 meters of the tube. This measurement agreed to 
within 1 percent with the value obtained from pressure and phototube measurements along 
the tube. 


Density Parameter 

The absorption technique yields a measurement of the integrated density of O 2 along 
the absorption path. This O 2 density can be related to a total density but this relationship 
requires the use of assumptions about the structure of the test ’’slug.” Rather than impos- 
ing these restrictions the density will be expressed in the parameter p^. The total den- 
sity is then obtained by using 


P = 



(7) 


This equation requires the choice of an l that corresponds to an effective flow diameter; 
and the mass fraction F of O 2 in the flow. The range of proper values for l and F 
will become apparent when the data trends are examined. 

It is now possible to estimate the pitot pressure based on this density and the meas- 
ured interface velocity. Thus, 


t,c 


P P 
° 2,0 d 
DF 


V 


( 8 ) 
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For these initial conditions the pitot pressure drops sharply, or dips, during the test 
flow. One objective of this study was to obtain more complete data for the dip condition 
to aid in explaining this phenomenon. 

Measurements during this dip condition are shown in figure 9; initial conditions are 
listed in table 1. Pitot pressure p^ starts well below the ISFE prediction and increases 
for about 100 microseconds. (See fig. 9(a).) During the period between 40 and 165 micro- 
seconds, p^ lies between the ISEE and ISFE calculations. A sudden drop in p^ ends 
this period of "uniform" flow. After the dip occurs, the pressure record becomes very 
spiked but it is clear that the dip lasts about 100 microseconds. Scale depth p^ (fig. 9(b)) 
shows the same general trend as p^.. Values first increase to a reasonably steady level 
and then decrease. This decrease starts earlier and is slower than that for p^. This 
difference is thought to be caused by the integrating effects of the p^ measurement. It 
is interesting to note that m the p t dip, p^ is reduced by about 50 percent, roughly the 
same percentage as occurs in p^. This agreement implies that the pitot pressure dip is 
a direct result of density decrease in the flow. Values of T y (fig. 9(c)) start near 
the values of the ISFE calculation, fluctuate a bit, and then decrease slowly for about 
40 microseconds. The temperature i]emains uniform with only minor, short-lived devi- 
ations for the rest of the test time. These deviations probably result from the data- 
reduction technique rather than from real changes in T y q . Changes in A (fig. 9(d)) 

are noticeably different for this example. There is a general increase, lasting 40 micro- 
seconds, followed by a uniform level at 1.5 ± 8 percent, much larger than the 0.8 value in 
the first example. The sudden changes that occur at the time of the dip probably indicate 
a breakdown in the assumption of uniformity in the radial direction. 

Data Trends 

As a means of obtaining data trends, tests have been grouped according to test gas 
composition, p^ and Vg j. The results of each group were then plotted as a function 
of V^. Since the primary diaphragm was pressure burst, the scatter in driver pressure 
caused a small random variation in Vg j. In order to correct for this variation, all data 
for the same nominal conditions behind the primary shock were standardized to an averagf 
Vg i and the data were corrected accordingly. Thus, a measured value Y (pressure 
or temperature) was corrected by using the relation 

/ Y ISEE 2 Y corr\ 

\ y isee ' y efe/ 

V s,s Y s,m 


Y ISEE ~ Y meas 
Y EEE " Y ISFE 
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78 percent N 2 mixture at the initial conditions of example 1 are shown in figure 13. 

These data were generated from a computer program originated by C. J. Schexnayder 
and J. S. Evans of the Langley Research Center. This program solves for the flow 
behind a normal shock. It allows chemical nonequilibrium but does assume vibrational 
equilibrium. In figure 13 the abscissa is the time of entrance of the particle path into 
the expansion fan and the ordinate is the ratio E of chemical energy tied up in disso- 
ciation to the corresponding value for equilibrium conditions. Thus, if equilibrium is to 
be attained, this ratio must be unity. That part of the flow that will form the test slug 
enters the fan during the first 30 microseconds. Consequently, according to this calcula- 
tion the flow in region 2 should be considered to be out of equilibrium, and nonequilibrium - 
flow calculations for this type of expansion tube should start at the primary shock - not 
at the secondary interface as in references 6 , 7, and 8 . 

It should be noted that, because of the nature of the expansion-tube flow cycle, a 
significantly lower level of dissociative nonequilibrium is achieved than in a steady flow 
nozzle for the same test condition. The expansion tube avoids the high enthalpy state of 
a stagnation reservoir, so that the peak temperatures produced in such a reservoir are 
reduced. Some perspective of the present level of dissociation is obtained by comparing 
the highest value of the compressibility factor Z in the present study (corresponding 
to equilibrium behind the primary shock), which is less than 1.04, with the equilibrium 
value of 1.45 which would occur behind a bow shock on a model for this test condition. 

A summary of A is shown in figure 14; A was defined as 


A = 


O 


2,0 




DFp 


(ID 


Several assumptions are included in this definition and the breakdown of these assumptions 
will cause this parameter to be different from unity. For example F has been assumed 
to be F 2 , but according to the results detailed in figure 13, F should be closer to F^ 
(about 15 percent larger for the average run). Variations of V 5 with time are also 
important, but the largest error is probably caused by a nonuniform radial distribution 
of temperature, density, and flow velocity. 

The method of illustrating trends in A in figure 14 is different from that used for 
other variables. The early value is still indicated by the base of the arrow, but the head 
of the arrow is replaced by the symbol and gives the value of A 75 microseconds back 
in the flow. In all cases A had reached a steady value (±10 percent) by this time and 
this value represents the test core value. 

It may help to remember the manner in which the test conditions vary with interface 
velocity; that is, for the same conditions behind the primary shock, expansions to higher 
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The wall pressure and pitot pressure measurements, in accord with those obtained 
in previous studies in this facility, approximated the value calculated by the incident- 
shock, frozen- expansion model or fell between this value and the incident- shock, 
equilibrium -expansion calculation. None of the measurements suggested that any 
improvement in the calculation flow model was achieved by including the effects of a 
reflected or standing shock at the secondary diaphragm. The wall pressure, pitot pres- 
sure, and temperature measurements indicate that the simplest reasonably accurate pre- 
diction of test gas state would neglect any shock reflection at the secondary diaphragm 
and use the mean of the equilibrium and frozen expansion models. 

Density measurements at conditions for which the pitot-pressure dip occurs indicate 
that the pitot pressure decrease is principally the result of a density decrease. 

The path integrated density measurements at conditions free from the pitot-pressure 
dip are, of course, influenced by the wall boundary layer. The results of calculations 
show that these densities can be brought into agreement with the measured pitot pressure 
and interface velocity by assuming that the boundary layer accounts for a 20-percent 
decrease in the effective path length. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 12, 1973. 
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APPENDIX A - Concluded 


The error in a from the relation 


a 


A 147 

A 160 


becomes 


A a 



1/2 



This Aa was used with the data of figure 4 to calculate the error involved in 
measuring T y by assuming a 1-percent error in the measurements of 1^ Q and 
1^ for each wavelength. The results of these calculations are shown in figure 16. The 
data covered in this study fall in the range 2 < < 12. 
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Expansion-interface reaction 



Expansion-tube flow cycle showing pertinent flow regions 







Wall pressure 


i 


Flow 


Trigger 


High 

voltage 

trigger 


— >_ pq— 4A 




U 
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Figure 9.- Example of typical test results at the pitot pressure dip 
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Figure 12.- Continued 
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